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1. Initialize the loop 



Ipl = { 



{LCLi, LCLj) , {LCLir UCLj) , 

{UCLi, UCLj) , {UCLi, LCLj), (LCLi, LCLj) 



k - 



1. 



2. Create the half space hs defined by the specification k: 

hs = {X G I fi,(X) > LSLj,] . 

For two-side specifications, repeat Step 3^5 for USLk. 

3. Set IpO = Ipl, 

~ vi, and 

Vjb = V2. 

Empty Ipl; 

4. If both and Vj^ are inside hs then 

add vt into Ipl; 
else 

if one of Va and vt, is inside hs then 

calculate the intersection point Vc of VaVjb and hs 
add Vc into the new loop ipl 

if is inside hs then 



5. Set Va - Va+i and 

If not all vertices are done, go to step 4. 

6. A- = A- + 1. If A" < n, go to step 2. 

7. Add the first vertex into the end of Ipl to finish the 
loop. The region surrounded by Ipl is the feasible space 
of Xi vs. Xj with other parameters fixed. 



add into Ipl . 
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1. Choose n distinctive constraints from all m 
specifications and n parameter limits; 

2. LU decompose F into L-U; 

3. Constitute one distinctive right-side Y from, 
specifications and parameter limits. 

4. Solve LUX = Y by forward and back substitution. If the 
solution X* satisfies all the specifications, add it into 
the extreme point list. 

5. If all distinct Ys have been operated upon, go to next 
step. Otherwise go to step 3. 

6. If all distinct constraint combinations have been 
operated upon, go to next step. Otherwise go to step 1. 

7. All extreme points and the facets are now available. The 
feasible decision space and performance space can be 
acquired by projecting the facets into each 2-dimensional 
viewing pane, though other projections and 
representations may be beneficial. 

8. The boundary of the feasible space in each viewing pane 
is a convex hull of the extreme points. Graham-scan 
algorithm can be adopted to retrieve the convex boundary 
in 0{hlq{b)) time where h is the number of all vertices. 
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